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We propose a protein model based on a hierarchy of constraints that force the protein to follow 
certain pathways when changing conformation. The model exhibits a first order phase transition, 
cooperativity and is exactly solvable. It also shows an unexpected symmetry between folding and 
unfolding pathways as suggested by a recent experiment. 
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Proteins function by switching between different well-defined conformations in the native state. In order to reach 
' the native state the protein has to fold. The number of possible conformations is huge, and the existence of folding 
^\ I pathways guiding the proteins between the relevant conformations seems necessary The existence of pathways 
' has been experimentally established [^|-^, but the mechanisms behind these are still to be understood. The free 
(— I , energy difference stabilizing the folded state of a protein, is of the order 20 kBT^oom for a typical domain consisting 
a ■ of about 100 amino acids, meaning that the binding energy per amino acid is only of order of kBTroom §• This is a 
' ' , surprisingly small value. The mechanism responsible for stability at very little cost is not known, but has been named 
''si" ' "cooperativity" ||]. 

' , In this letter we propose to model protein folding by a parametrization of the folding pathway. Cooperativity 
' appears in the theory through assuming that only interlocked degrees of freedom are active in the energy window that 
J> . separate the folded and unfolded states. This also leads to a first order folding transition. 

' Furthermore, the model exhibits an interesting dynamical symmetry when changing conformation. One would 
t , expect it to follow the folding pathway in opposite direction when it unfolds. However, under certain conditions it 
moves in the same direction along the folding pathway both when folding and unfolding. We call this the "mirror 
effect." This feature of the model is in agreement with recent measurements |^ which suggest this mirror symmetry 
' in the conformational change of albumin at low pH. 
On ■ The amino acids constituting a protein may rotate relative to each other, thus inducing conformational changes in 
the protein. They furthermore interact. Through changing bond angles, amino acids far apart on the protein may be 
brought close to each other to form bonds that stabilize a given conformation. 

Simple models for protein folding which encompass guiding |Q and cooperativity |^ have been proposed. More 
^ complex models based on spin glass hamiltonians |^ need to introduce "folding funnels" in the energy landscape in 
^ order to guide the slow dynamics of a frustrated system . The existence of a pathway implies that the bond 

Q , angles change in a certain order, i.e. there is a hierarchy among the bond angles themselves. In order to clarify the 
O ■ significance of such a hierarchy, we introduce a model that simultaneously displays a folding pathway, cooperativity 
and a first order phase transition, thus reproducing three important characteristics of real proteins. The model may 
be visualized as in Fig. ||. The variables are effective folding angles ipi to ipiq, each associated with matching a given 
sub part of the protein into its position in the ground state (the template). The key point of the model is that a 
■ correct fold at level i gains binding energy when all lower levels are correctly folded. Our angle variables are not to be 
thought of as the (f), ip angles of the polypeptide backbone, but rather as coarse grained variables (i.e. averaged over 
a small number of monomers). An intuitive picture would be the angle between two a helices connected by a loop. 
Our model corresponds to the interaction Hamiltonian 

H = -ai ifii - 02 ipiip2 - as (pi(p2V3 cln ipi(p2 ---^n , (1) 

where ai to a^r are coupling constants. 

In order to present explicit calculations we scale all angle variables ipi to be between and 1 and set all coefficients 
Ui equal to 1. Furthermore, we assume for the moment that angle variables only take discrete values, and 1. The 
Hamiltonian then represents the number of folded variables {tpi — 1) counted from i = 1 until the first i = rif + 1 
where ifi = 0, and can simply be written 
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H = -nf 
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To explore the thermodynamic properties of the model, we calculate its partition function by summing over degen- 
erate states, 

Z{T,N) = y 2^-i-"/e"//^ + e^/^ = 2^-i exp(jV(l/r - In 2)) - 1 ^ ^^^/t ^ (3) 
^-^ exp(l/r — ln2) — 1 

where T is temperature measured in units of the coefficients ai — 1. In the limit N 00, the partition function has 
a discontinuous derivative for T = Tc = 1/ In 2. This means that the system undergoes a first order phase transition 
at this temperature. At temperatures T < Tc the binding energy of the folds dominates and the system freezes into 
the folded state with (n/) f» N, whereas for T > Tc the winning term will be the entropy that is gained from each 
folded variable, yielding (n/) ~ 0, corresponding to an unfolded, molten state. Experimentally, it is known that small 
proteins typically undergo a first order transition, while for large, multidomain proteins the situation can be more 
complex [||. 

In addition to a phase transition, the model exhibits cooperativity, manifested through the fact that below Tc, say 
for T = 1, the chain is folded even though the binding energy per link is only fcT, which would lead to a disordered 
(unfolded) state in the absence of cooperativity. In the present model, at T < Tc, nearly all states will match the 
template, with fluctuations only at the bottom of the hierarchy. This matches the findings from real proteins, where 
cooperativity implies that most degrees of freedom are locked by others. 

During a slow "adiabatic" change of temperature through the phase transition, the folding/unfolding follow specific 
pathways. Under cooling, the first level in the hierarchy must fold first, 991 = ^ 1 in order for the energy to change 
at all. This is due to the nature of the Hamiltonian where ipi is present in all the terms, see (|l|). After the first, 
the next level in the hierarchy, ip2, must fold in order to change the energy further and so on. Under heating, the 
unfolding will take place from bottom and up, i.e. first ipN = 1^0 leading finally to ipi = I —^ 0. This is due to the 
fact that the smallest increase of energy is given by the unfolding of the deepest levels in the hierarchy. 

Under non-adiabatic changes the unfolding is however different. When T < Tc, the free energy is 

Fe{T,N) = ~N, (4) 

as it is dominated by the binding energy. For T > Tc, the free energy is 

FsiT,N) ^ ~TNlii2 , (5) 

which results from entropy now dominating. We denote Fe the energy branch and Fs the entropy branch of the free 
energy. We may envision these extended beyond the critical temperature Tc as shown in Fig. ^. We now perform a 
quench of the system, i.e., rapidly change its temperature. Starting in equilibrium at a temperature T < Tc, marked 
"1" in Fig. H, and rapidly heat the system to T > Tc, it moves to "2" in the figure, following the energy branch of the 
free energy. The system then thermalizes by moving to "3" , after which the system ends up on the entropy branch. 
Starting with the system in equilibrium at a temperature T > Tc, point "4" in the figure, a rapid quench to a low 
temperature brings the system to "5" along the entropy branch. The system then moves from "5" to "6" to reach the 
energy branch which at this temperature is the lowest free energy. 

The question is now how the angle variables react to these two types of quench. When the system is rapidly cooled, 
the angle variables have to fold in the order 1,2,3,... N — otherwise, no change in energy occurs. This is the same 
order as the system follows during adiabatic cooling. However, when going from low to high temperature the situation 
is different. When the system moves from "2" to "3" in Fig. |^, the unfolding rate, F^, for each ipi, is dominated by 
the size of the phase space of the final state, 

Therefore the rate increases with decreasing i, i.e., the rate of change is largest for (pi, second largest for (p2 etc. This 
is the same order as when the system was rapidly cooled and the opposite to the order found in adiabatic melting. 
We call this reversal of direction "mirror effect." So far as we have regarded the variables ipi as binary, as soon as 
the first angle has unfolded, the others will unfold in no particular order, see Fig. |l]. With continuous variables this 
degeneracy of the model is lifted. 

In a recent experiment, an indication of the existence of this "mirror effect" was found In the experiment, 
the conformational change induced in albumin by lowering pH is studied on samples consisting of very few (1-100) 
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molecules. The technique is described in detail in briefly, the molecules are sandwiched between a flat plate and 
a micron size sphere, and the movement of the sphere relative to the plate, measured to sub nm resolution with 
an evanescent wave technique, is used to reveal changes in shape (e.g. swelling or contracting) of the molecules in 
response to a fast pH change. This corresponds to quenching the proteins into the new state. 

Albumin is a large protein (~ 600 aminoacids), consisting of 6 domains. It is believed that the conformational 
change at low pH consists essentially in a rearrangement of these domains, so that the shape of the molecule changes 
(becoming about 60 % longer). However, this is not a complete unfolding: the molecule is still compact, and elements 
of secondary structure persist. In Fig. ^ we present some measurements from the experiment Q. As described 
in 1^, the conformational change is seen to proceed through a series of ^ 5 steps in time, suggesting a sequential 
rearrangement of the domains. Thus for this process, a well defined pathway exists. 

The experiment further suggests that the reverse conformational change (obtained by restoring neutral pH con- 
ditions) proceeds through the same sequence (not the time reversed sequence) of steps. We observe that the steps 
occur at similar positions, supporting a mirror imaging of the retracing pathway. In support of this observation, we 
project in Fig. ^ the trajectory on the distance axis, to show a histogram of the time spend at each distance value, 
i.e. such that flat regions (waiting times) will be represented by peaks in the histogram. We can then average several 
histograms. The comparison of the two symmetries, i.e. the "mirror effect" (b) and simple time reversal (c), supports 
that the proteins under these conditions show the "mirror effect." 

In conclusion, we believe to have introduced the first model which simultaneously reproduces three key features 
of real protein folding, namely a folding pathway, cooperativity and a first order folding transition. In addition the 
model raises some interesting symmetry discussions (i.e. the "mirror effect"). 

A.H. thanks Nordita for funding his stay in Copenhagen during which this work was done. Also we like to thank 
R. Donangelo and K. B. Lauritsen for discussions. 
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FIG. 1. The figure describes the refolding of a connected structure, such as a partly unfolded protein, into a well defined 
template. It illustrates schematically how one by successive foldings of angle variables can match the edges 1,2,3.. to the 
template edges l',2',3'... We suggest to copy the figure and do the folding. One will see that it is not possible to match say 
edge 2 to template 2', without first matching 1 to 1'. In general to match j with j' demands correct matchings for all i < j. 
This defines a unique pathway from the unfolded to the folded structure, in a way that can be quantified by the Hamiltonian 
in Eq. 



FIG. 2. The free energy per variable. The full line show the equilibrium free energy as a function of temperature. The 
dashed lines show the entropic and energetic branches of the free energy, respectively (Eqs. (^) and (^)). Path 1^2 shows a 
rapid heating of the system. It then unfolds following path 2 — > 3. Path 4 — > 5 shows a rapid cooling of the system followed by 
the folding following path 5 — » 6. 
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FIG. 3. Curves from the experiment ||6| showing steps in the dynamics of the forward (left curve) and backward (right curve) 
conformational change. The ordinate is the height of the sphere above the plate, h, in nm, and the abscissa is the time in 
seconds with arbitrary zero point. For the right curve, only the distance axis has been reversed. The two curves have been 
placed next to each other for ease of comparison. The dashed lines help identify the steps. The correspondence between the 
steps in the two curves is the "mirror effect." 

FIG. 4. Histograms of the time spent at each distance value, constructed from projecting curves of the type shown in 
Fig. ^ on the distance axis. The waiting time regions in the original curves then give rise to peaks in the histogram and the 
distance between two peaks is a step size. We then average several histograms, (a) is an average over histograms corresponding 
to three unfolding curves. The large peaks at 8 and 11 nm correspond to the two baselines in Fig. pi The four peaks in 
between correspond to the four waiting times shown by the dotted lines on Fig. ^, left curve, (b) is an average over histograms 
corresponding to three folding curves transformed by reversing the distance axis (i.e. according to the "mirror effect"). Here 
there are only three peaks corresponding to the three steps in Fig. ^ right curve. As can be seen by comparing histograms (a) 
and (b) the missing step is either covered by the peak of the base line or simply does not occur in refolding, (c) presents the 
same histogram as in (b) but with simple time reversal, i.e. no reversal of the diameter axis. We see that (b) leads to better 
overlap with (a) than (c), indicating that the proteins under these conditions follow the "mirror effect." 
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